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p ^ , Abstract 

I We consider Schrodinger equations with variable coefBcients and 

• ■ the harmonic potential. We suppose the perturbation is short-range 

I type in the sense of [7]. We characterize the wave front set of the 

■ solutions to the equation in terms of the classical scattering data and 

the propagator of the unperturbed harmonic oscillator. In particular, 
we give a "recurrence of singularities" type theorem for the propagation 
of the period t — tt. 

> 

; 1 Introduction 



In this paper we consider a Schrodinger operator with variable coefficients 
and the harmonic potential: 



in 



00 , 1 '" \ 

O ■ ^ = ~2 ^ d^,ajkix)do,, + + V{x) 

^ . on IK = L (M"), n > 1. We denote the unperturbed harmonic oscillator by 

!=! : Ho: 

Hq = --A + -\x\^ on^K, 
and we suppose H is a short-range perturbation of Hq in the following sense: 

Assumption A. ajk{x), V{x) G C°° (M" ; M) for j, A; = l,...,n,and {ajk{x))j^k 
is positive symmetric for each x E M". Moreover, there exists /i > 1 such 
that for any a G 



|a^(a,-fc(x)-5,fe)| <C,(x)-^-H, 
\d'iV{x)\ < C«(x)2-^-l"l 
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for a; G M" with some Cq > 0, where (x) = ^1 + 

Then it is well-known that H is essentially self-adjoint on Cq°{MP), and 
we denote the unique self-adjoint extension by the same symbol H. We 
denote the symbols of H, Hq, the kinetic energy and the free Schrodinger 
operator by p, po, k and ko, respectively. Namely, we denote 

1 1 11 

Pix^0 = 2 ajk{x)^j^k + ^\xf + V{x), Po{x,0 = ^\^\'^ + ^\x\^, 

j,k=i 

1 " 1 
Hx,0 = 2 ^ (^jk{x)Cj^k, ko{x,^) = 
j,k=i 

We denote the Hamilton flow generated by a symbol a(x,^) on M^" by 
ex.p{tHa) : M?"" M^". We also denote 

7ri(X) = X, Tr2{X) =^ for X = {x, ^ e 

Let (xo,^o) G I^^"- {xo,Co) is called forward (backward, resp.) nontrap- 
ping (with respect to k) if 

|7ri(exp(ti7fe)(xo,Co)| ^ oo 

as t — > -|-oo {t — > — oo, resp.). If (xo,^o) is forward/backward nontrapping, 
then it is well-known 

{x±,C±) =1™ exp(-ti?fc ) oexp{tHk){xo,Co) 

t — >ihoo 

exists, and S±: (xo,^o) '-^ {x±,$,±) are locally diffeomorphic (see, e.g., Naka- 
mura [7], Section 2.). 

Now we present our main results of this paper. Let us recall our harmonic 
oscillator Ho has a period 27r, i.e., e~^'^'^^°(p = ip fov (p E Ji. Moreover, we 
have 

e^*^-^V(x) = v9(-x), ip£j{. 
Our first result concern the evolution by H up to time tt. We denote 

u{t) = e~^^^uo, uo e IK. 

We denote the wave front set of a distribution / by WF{f). 

Theorem 1.1. (i) Suppose (xo,Co) is backward nontrapping, and let < 
to < TT, Uo & IK. Then 

(xo,Co) e WF{u{to)) ^ (x_,C-) e Wy(e-^*°^°no). 

(ii) Suppose (xo,^o) is forward nontrapping, and let — tt < to < 0, uq £ Oi. 
Then 

(xo,eo) G WF{u{to)) ^ (x+,e+) G W(e-'*°^°uo). 
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Remark 1.2. We note that microlocally e ° is a rotation in the phase 
space. More precisely, for any reasonable symbol a = a{x,^), 

e-^*^°a"'(x, D^)e'^"° = a'"{cos{t)x + sm{t)D-„ - sm{t)x + cos{t)D^), 

where a''"(x,Dx) denotes the Weyl-quantization of a. Hence, in particular, 
(xo,Co) e WF{e-'*"°uo) if and only if there exists a symbol: a G C^(R2") 
such that a{xo,^o) ^ and 

\\a'"{cos{t)x - sin(t)i:>a;,/i(sin(i)a; + cos(t)i:'^))uo|| = 0{h°°) 

as h ^ 0. 

At the time t = ivr, u(t) behaves differently. We denote the set of 
forward (backward, resp.) nontrapping points by T+ (T_, resp.) C M^'^\0 := 
{{x, e C 0}. S± are diffeomorphism from (R^" \ 0) \ T± to R^" \ o, 
and hence S^^ arc well-defined from R^"\0 to (R^"\0)\T±. We also denote 
the antipodal map in R^" by F, i.e., r(x,^) = (—a;, — ^). 

Theorem 1.3. (i) Suppose (xo,^o) is backward nontrapping, and letug G "K. 
Then 

(xo, Co) G WF{u{tt)) ^ S-^oTo S_(xo, Co) G T^F(«o). 
('nj Suppose {xo,^o) is forward nontrapping, and let uq G "K. Then 
(xo,eo) G WK-tt)) ^ 5zioro5+(xo,eo) G 

The microlocal singularities of solutions to Schrodinger equations have 
been attracted attention during the past years, especially after the publica- 
tion of the break-through paper by Craig-Kappeler-Strauss [1] in 1996 (for 
more literature, see references of [1] and [7]). On the other hand, the sin- 
gularities of solutions to the harmonic oscillator type Schrodinger equations 
have been studied by several authors, including Zelditch [11], Yajima [9], 
Kapitanski-Rodnianski-Yajima [6], Doi [2] and Wunsch [10]. Most of these 
works concern the case with constant coefficients with potential perturba- 
tions. In particular, if the metric is flat, i.e., if ajk{x) = Sjk, then 5± is the 
identity map, and Theorem 1.3 recovers results in [11], [9], [6] and [2]. In 
fact, if the metric is flat. Theorem 1.1 is also obtained by Doi [2], and this 
paper is partially motivated by this beautiful work. Moreover, he also con- 
sidered a class of long range type perturbations, i.e, when V{x) = 0{\x\) as 
|x| — >■ oo, and demonstrated that a shift of singularities occurs. Microlocal 
smoothing effect for the Schrodinger equations on scattering manifolds with 
harmonic potential is studied by Wunsch [10]. 

In a sense, this paper is an analogue of a work by Nakamura [7] where the 
microlocal singularities of asymptotically flat Schrodinger equations is stud- 
ied (see also a closely related work by Hassel- Wunsch [3] and Ito-Nakamura 
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[5]). The main difference (and the novelty) is the analysis of the classical tra- 
jectories with high energies. In [7], the standard classical scattering theory is 
sufficient to prove the propagation of singularities (with a scaling argument). 
However, in the presence of the harmonic potential, the high energy asymp- 
totics of the classical trajectories is completely different from the long-time 
asymptotics. Thus we need to obtain precise high energy asymptotics of 
the trajectories using the time evolution of the harmonic oscillator, and it 
is carried out in Section 2. This situation is somewhat similar to the case of 
such analysis for Schrodinger equations with long-range perturbations [8], 
but the asymptotics itself is naturally completely different. We also note 
that our results have much in common with a paper by Zelditch [11], at 
least in spirit, and our results may be considered as generalizations of his 
results to variable coefficients cases. 

The main results are proved in Section 3, and the argument is similar to 
[7]. However, the scaling argument is slightly more complicated, and we try 
to give a more transparent formulation (see the last part of Section 2 and the 
beginning of Section 3). In the last section, we discuss generalizations of our 
main theorems to the case when the harmonic potential is inhomogeneous. 

Throughout this paper, we use the following notations: Our Hilbert 
space is "H = L^(M"), and the space of the bounded operators on J{ is 
denoted by £-(CK). The Fourier transform is denoted by 

uiO = 9^u(0 = (27r)-"/2 f e-'-<uix)dx, 

and the inverse Fourier transform is denoted by u(x) = 3"*u(x). For a 
smooth symbol a{x, ^) on R^", we denote the Weyl quantization by a^{x, D^), 
i.e., for u G S(M"), a Schwartz function on R", 

a-(x,r»,)u(x) = (27r)-"|y" e'(--yy^ai^,Ouiy)dy dC 
We use the S{m, g)-symbol class notation of Hormander [4]. 

2 High energy asymptotics of the classical flow 

In this section, we study the high energy behavior of the classical flow gener- 
ated hy p{x, ^). More precisely, we consider the properties of exp(ii/p)(xo, A^o) 
as A ^ -l-oo. Throughout this section, we suppose (a^o^^o) is forward non- 
trapping, and consider the case t > 0. The case t < can be considered 
similarly. 
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For A > 0, we write 

1 " |xP 1 

j,fc=i 

Then, by direct computations, we learn 

(2.1) 7ri(exp(iiJp)(x, XO) = 7ri(exp(AiiJpA)(x, 0), 

(2.2) 7r2(exp(iiJp)(x, A^)) = A • 7r2(exp(At/fj,A)(x,e)). 

Hence, it suffices to consider e^p{tHp\){x,$,) for < t < Xto, instead of 

exp{tHp){x, XC) for < t < to- 
We note, for each fixed t G M, 

(2.3) lim eMtHpx)ix, ^ = exp{tHk){x, ^ 

by the continuity of the solutions to ODEs with respect to the coefficients. 

Hence, if t > is large and then A > is taken sufficiently large (after 
fixing t), 7ri{cxp(tHp\)(xo, $,o)) is far away from the origin by virtue of the 
nontrapping condition. The next lemma claims that this statement holds 
for < t < A(5 with sufficiently small 6 > 0. 

Lemma 2.1. There exists 6 > and a small neighborhood J7 o/(xo,Co) such 
that 

|7ri(exp(tiJpA)(x,0)| > cit -C2 forO<t< XS, {x,0 G ^ 
with some ci,C2 > 0. 
Proof. In the following, we denote 

exp{tHpx){x,0 = {y^{t;x,C),7]^{t;x,0)- 

By the conservation of the energy: p'^{y'^{t),r]^{t)) = const., and the ellip- 
ticity of the principal symbol, we easily see 

^\y\t;x,Of + \v\t;x,0\^ < C, {x,0 en,teR, 

where is a small neighborhood of (xo,^o)- Hence, in particular, we have 

\y\t;x,0\<C{t), \ri\t;x,0\<C 

for t > with some C > 0. On the other hand, by direct computations, we 
have 

= 4p\y\ri^) + 2W{y\r]^), 
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where 



j,k,e,m 



j,k,e,m j,k 



Combining these, we learn 

We note p'*'(xo,^o) = k{xo,^o) + 0{\~^) and k{xo,^o) > 0, and hence 
pHxo,^o) > for large A. Since \-'^{y^f = 0{{t/\)^), if < t < 5A 
with sufficiently small S > 0, the last term of the right hand side of the 
above inequality is small and 

^\y\t)\'>3p\y\v')-c,{y')-^ 

for the initial condition (x,^) G 0. By the nontrapping condition and (2.3), 
if To > is sufficiently large and A is large (depending on Tq), then 

C4(y\ro))-'^ < p\x,0 for (x,0 G and j^\y\To)\ > 0. 
Then by the standard convexity argument, we learn 

\y\t)\^ > \y\To)f+p\x,m - nf for t G [ro,^], 
and this implies the assertion. □ 

Lemma 2.2. Let 5 > Q and O as in the previous lemma, and let a G (0, 5). 

Then 

lim exp(— crAi? a) o exp((TAi?„A)(a;, ^) = Sj^{x,^), 
for {x,Oen. 
Proof. We denote 

iz\t; x, 0,C\t; X, 0) = expi^tHpx) o exp(ti7pA)(x, 0, 
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and we show the convergence of {z^{aX),(^{aX)) to S+{x,$,) for (x,^) G Q. 
We recall 

eM-tHpx){x,0 = (cos(l)x - Asin(|X, isin(i)x + cos(|)e) 
since Pq is the scaled harmonic oscillator. Thus 

z\t) = cos(|)y^(t) - Asin(i)r?^(t), = j.sm{{)y\t) + cos{{)ri\t). 
By direct computations, we have 



i4 = -\Mi)y'k+cos{{)^-cos{{)r,',-XsHj:)^ 

= -T sin(^)y^ + cos(^) ^ a,fe(2/^)??j^ - cos(^)77^ + \ sm{{)yj^ 



1 / A- 



= o{{t)-n 

for < i < 5A. Similarly, we have 



3 

i 

= 0(A-^(y^)-'^) + 0((y^)-''-^ + A-2(y^)-('^-i)) 
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for < t < SX. Moreover, for each t G M, we have 

3 «j 



^iVk-tfik), 



d 



where (y(f), 77(f)) = exp(fiJfe)(a;, ^). By using the dominated convergence 
theorem, we conclude 

hm z ((tA) = X + hm / —-—at = x+ I —{y — trjjat 

A— >oo A— >oo Jq at Jq dt 

= hm (y(f)-t?7(f))=7ri(5+(x,0), 

fcrA jaA /-OO J 

lim C^(ctA) = ^ + lim / ^dt = i+ -f,{t)dt 
= hrn fj{t) =7r2(S'+(a;,C)). 

T— >+00 

This completes the proof of the lemma. □ 

Lemma 2.3. Let < a < n, and let Q be a small neighborhood of (xo,^o) 
as in the previous lemmas. Then 

hm exp(— aAiJ a) o ex.p{aXH \){x, = S'+(x,^) 
A— ►oo ^0 

for {x,0 € n. 

Proof. It suffices to consider the case 5 < a < tt, and we fix such a. Let 
£ > 0, and we show that if 

ma^{\z\aX) - x+\, \C^{(7X) - C+|) > e 

then A is bounded from above, where (x+,^+) = 5+(x,^), and {z-'^X^) is 
as in the proof of the previous lemma. Our claim then follows from this 
assertion by contradiction. 
We first note 

y\t) = cos{{)z\t) + Xsm{{)C\t). 
For the moment, we suppose \z^{t) — x+\, \C^{t) — ^+| < e. Then we have 

\y\t)\ > Asin(i)(|^+| - e) - (|x+| +e) > S^X 
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with some 6i > provided SX < t < aX, e < |^+|/2 and A > Aq, where 6i 
and Ao depend only on and |^+|. Then, by using formulas (2.4) and 
(2.5), we learn 



dt 



z\t) 



< cx-^, 



dt 



Cit) 



< cx-^'-\ 



Now we choose A sufficiently large that 



maxi\z\6X)-x+\,\C\SX)-^+\)< 



2' 



and suppose 



max(|z^(aA) - x+|, \C\aX) - $+1) > s. 
Then there exists to € {SX, aX) such that 

msxi\z^{to) - x+l, \C\to) - e+l) = £ 

and 

max(|2^(t) - x+l, IC^(t) - ^+1) < £ for 5A < t < to.. 
By the above observation, we learn 



\z^{to)-x+\ = 



\SX)-x+ + j^^ ^{t)dt 



< \z^{8X) - x+l + C{to - 6X)X-f' 

< \z^{5X)-x+\ + C{a-S)X 



-(p-i) 



IC'(^o)-e^ 



c\6X)-c++ r 

Jsx 



to 



dt 



{t)dt 



<\Q\5X)-^+\ + C{a-5)X-^. 



Thus we have 



e = m^{\z\to) - x+l, IC^to) - C+l) <^ + C{a- S)X 



-(/x-l) 



and hence A < {2C{a — 5)/eY^^^ ^\ This proves the assertion. 
The next theorem follows immediately from Lemma 2.3. 



□ 



Theorem 2.4. (i) Suppose (xq, ^o) is forward nontrapping, and letO < a < 
TT. Then there exists a neighborhood Q of (xo,Co) such that 

lim 7ri(exp(-(7i7po) oexp((7iJp)(x,A0) = 7ri(S'+(x, 0), 

A— >oo 

lim A~V2(exp(-crifpo) o exp(crifp)(x, A^)) = 7r2(S'+(x, ^)) 
A— >oo 
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for (x,^) G fl, and the convergence is uniform in fi. 

(ii) Suppose {xq,^q) be backward nontrapping, and let — tt < a < 0. Then 
there exists a neighborhood o/(a;o,^o) such that 

lim 7ri(exp(-cri?po) oexp((7iJp)(x,A0) = 7ri(5_(a;, 0), 

A— >oo 

lim A"V2(exp(-cri7pJ oexp(cri7p)(x,A0) = 7r2(S'_(x, ^)) 

A — *oo 

for {x,$,) G il, and the convergence is uniform in CI. 

We introduce several notations as a preparation for the proof of our main 
results. We set 

l{t-x,$,) = {p o exp{tHp^)){x,0 - po{x,0 

n 

= (ajfc(cos(i)x + sin(i)^) — 6jk) (— sm(t)xj + cos(t)^j) x 
j,k=i 

X (- sm{t)xk + cos(t)^fc) + V{cos{t)x + sin(t)^). 
Then it is easy to show that i{t; x, £,) generates the scattering time evolution: 

St = exp(-tiJpJ o exp(iiJp). 

Similarly, 

e\t; X, = {p^ o eMtHpx)){x, - P^{x, 
generates the time evolution: 

= exp{-tHpx) o eMtHpx). 

We denote the scaling with respect to ^ hy 3\, i-e., 

dx{x,0 = {x,XO for(x,e)GM2" 

Then by (2.1) and (2.2), we have 

exp(tiJp) o3x = dxo exp(AiiJpA), exp(tiJpJ o 3^ = 3x0 exp(Ati?pA), 

and hence we also have 

(2.6) Sto3x = 3xoS^t. 

3 Proof of main results 

In this section, we mainly concern the case {xq, ^q) is forward nontrapping, 
and prove the part (ii) of Theorem 1.1. 
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We first consider the property of e^tHo^-itH f^j. ^ < q. Let vq G C^(R"), 
and we consider 

v{t) = e'^^'^e-'^^VQ. 

Tfien it is easy to observe 

= -i[e''"oH^-itHo _ Ho)v{t) = -iL{t)v{t), 

where L{t) = e«*^o/fg-»t^fo _ jj^ \Ye recall, for any reasonable symbol 
a{x, we have 

e^*^°a"'(x,L»^)e-^*^° = (a o exp(ti7pjr(x, 

without the remainder terms, since Po{x, ^) is a quadratic form in {x, and 
we employ the Weyl quantization. Thus we have 

L{t) = {pocMtH.p,)r{x,D^)-po{x,D,)=r{t;x,D^), 

i.e., i{t:x,^) is the Wcyl-symbol of L(t). This is in fact expected, since 
^itHo^-itH -g ^j^g quantization of St- 

Let O be a small neighborhood of {xq, Co) as in the last section, and let 
/ G C^(n) be such that f{xo, Co) > 0, and f{x, > on R^n ^e then set 

fh{x,0 = fix,hO, h = x-\ 

where /i > is our semiclassical parameter. We consider the behavior of 

G{t) = e^*"°e-'*"r{x,hD^)e^*"e-'^^° 

as /i — > 0. The operator valued function G{t) satisfies the Heisenberg equa- 
tion: 

(3.1) j^Git) = -i[L{t), Git)], G(0) = r (x, hD,). 

The corresponding canonical equation of the classical mechanics is 

^{t;x,0 = -{£,M{t;x,0, mx,0 = fixM) 

and the solution is given by 

Mt;x,0 = ifhoSi')ix,o 

since St is the Hamilton flow generated hy £{t;x,^). Now we note 
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Hence, recalling (2.6), we learn 

fh o Si' = / o o = fo{Sto 3x)-' 

In other words, we have 

V'o(«;x,0 = {fhosr'){x,0 = {fo{s^,)-'){x,hO. 

We expect 

G{t) ^ Mf, X, D,) = if o (S^tr'Kx, hD,) 

for small > 0, and we construct the asymptotic solution to the Heisenberg 
equation (3.1) with the principal symbol ^o(i;a^5 0- 

Lemma 3.1. Let — tt < to < 0, and set I = [to^O]. There exists ijj{t;x,^) € 
C^(R2») fort el such that 

(i) ^{Q;x,i)= f{xM)- 

(a) il>{t;x,^) is supported in St o3x{^) = 3\ ° S^f{il). 
(Hi) For any a,/? G Z^J., there is Cap > such that 

|5^afv(t;x,0| < Caph\^\ tei,x,^e M". 

(iv) The principal symbol oftp is given by tpo, i.e., for any a,/3 G Z", there 
is Cap > such that 

\d^d^m;x,0-Mt;x,0)\ < Caph'+\^\, t e i,x,^ e R". 



(v) If we setGit) = ilj'"{t;x,D^), then 
d 



dt 



G{t)+i[L{t),G{t)] 



0(/l°°) 



as h ^ uniformly int&I. 



Proof. Given the classical mechanical construction above, the construction 
of the asymptotic solution is quite similar to (or slightly simpler than) the 
proof of Lemma 4 of [7]. Wc note e{t; x, C) G S{{0'^,dx^ + d^^/{^f ) locally 
in X, and 7ri(supp (V'o(^j ') "))) contained in a compact set by virtue of the 
asymptotic property: S^^ 5+ as A ^ oo. We omit the detail. □ 

Now the proof of Theorem 1.1 is almost the same as the proof of Theo- 
rem 1 of [7], and we simply refer the reader to the paper. 
Finally, we show Theorem 1.3 follows from Theorem 1.1. 
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Proof of Theorem 1.3. We prove the part (i) only. We note 

(3.2) WF{e-'^'^/^^^°u) = WF{u), 

(3.3) WF{e'^''/^^^°u) = WF{u) = r{WF{u)). 

In fact, e^*('^/^)^" is the Fourier transform. (See also Remark 1.2.) 

We set {x', f ) = 5*+^ o r o S-{xo,Co) so that 

(3.4) (x_,e-) = r(xV,e;), where ix'^,e+) = S+{x',0- 

By Theorem 1.1 (i) with to = 7i"/2, and u{Tr/2) as the initial condition, we 
have 

(xo,eo) e WF{u{7r)) ^ (x-,e-) e WF{e-''^'^/^^"°u{Tr/2)) 

^ {x-,^^)eWF{u{7r/2)). 

We have used (3.2) in the second step. On the other hand, by Theorem 1.1 
(ii) with to = — 7r/2, and u{'it/2) as the initial condition, and using (3.3), we 
also have 

(x', CO G WF(uo) ^ «, C;) G WF(e^(-/2)^o„(7r/2)) 

= WF{u{tt/2)) = V{WF{u{it/2))). 

By (3.4), this implies the claim of Theorem 1.3 (i). The part (ii) is proved 
similarly. □ 

4 Inhomogeneous harmonic oscillators 

Here we consider the case when the harmonic potential is inhomogeneous, 
i.e., 

1 1 " 
Hq = --^ + - ^ bijXiXj, 

with a positive symmetric matrix (bij), and 

^ n \ ^ 

H = --Y^ dxiaij{x)d:cj + 2 X] ^iJ^i^j + ^i^)- 

i,j=l i,j=l 

We assume {ajk{x)) and V{x) satisfy Assumption A. By an orthogonal trans- 
form, we can diagonalize the harmonic potential, and hence we may assume 

^bijXiXj = Yl]=i^j^'ji where > 0, j = l,...,n, are eigenvalues of 
(bij). The behavior of the inhomogeneous harmonic oscillator depends on 
the number theoretical properties of (vj)^^^. If there exist no > such 
that 

(4.1) tQi^j G vrZ, j = 1,... ,n, 
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then it is well-known that the recurrence of the evolution operator does not 
occur, i.e., there are no to 7^ such that e~'*^° = /. In this case we have 
the following result: 

Theorem 4.1. Suppose (a;o,Co) is backward nontrapping, and suppose that 
there are no to > such that (4-1) hold. Then for any t > 0, 

ixo,^o) eWF{e-'^"uo) ^ {x-,^-) eWF{e-'^^°uo). 

Obviously, an analogous result holds for t < 0, but we omit it here. 
If there exists to > such that (4.1) holds, then we have the following 
result: 

Theorem 4.2. Let to > be the smallest positive number satisfying (4-1), 
and let mj = tovj/ir G Z. We set 

T{xi, . . . ,Xn,^l, ■ ■ ■ ,^n) = {o-iXi, . . . , CrnX„, (Ti^i , . . . ,(T„^n) 



for (x,^) € M where aj = 1 if nij is even, and aj = —1 if rrij is odd. 
Suppose {xo,S,o) is backward nontrapping. Then for < t < to, 

and 

(xc^o) e M/F(e-'*"^no) ^ 5;^ o f o 5_(xo, ^o) e WF{uo). 

The proofs of these theorems are similar to Theorems 1.1 and 1.3, and 
we omit the detail. We only note the fact that 



exp 

for ue L^{R), j = l,...,n. 
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